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Abstract: We investigate the effects of hydrodynamic nonlocal response 
in hyperbolic metamaterials (HMMs), focusing on the experimentally 
realizable parameter regime where unit cells are much smaller than an 
optical wavelength but much larger than the wavelengths of the longitudinal 
pressure waves of the free-electron plasma in the metal constituents. We 
derive the nonlocal corrections to the effective material parameters analyti- 
cally, and illustrate the noticeable nonlocal effects on the dispersion curves 
numerically. As an application, we find that the focusing characteristics of 
a HMM lens in the local-response approximation and in the hydrodynamic 
Drude model can differ considerably. Interestingly, sometimes the nonlocal 
theory predicts significantly better focusing. Thus, to detect whether 
nonlocal response is at work in a hyperbolic metamaterial, we propose to 
measure the near-field distribution of a hyperbolic metamaterial lens. 
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1. Introduction 

Hyperbolic metamaterials (HMMs), also known as indefinite media, enjoy a great deal of 
attention owing to their unique hyperbolic dispersion relations im ff4l . with associated high- 
wavenumber propagating waves without upper limit. This leads to numerous applications, such 
as enhanced-light interactions iHJQ, sub wavelength imaging jS HTSl . negative refraction [I], 
and low-loss fiber cladding llT4l . The HMMs are usually artificially made by periodic dielectric- 
metal structures, such as a ID dielectric-metal Bragg grating. To describe the optical properties 
of the HMMs, it is common to employ the local -response approximation (LRA), i.e. with every 
location r in the structure a certain value for the permittivity e(r) is associated. In the LRA, the 
effective material parameters of HMMs have been thoroughly studied ifTsHTgl . 

Thanks to advances in nanofabrication, we witness a miniaturization of the feature size of 
metamaterials towards the deep nanoscale. The LRA becomes more inaccurate, since the nonlo- 
cal response of free electrons starts to play a role [E0l - l32l . It is known that the nonlocal response 



causes a blueshift of the surface plasmon (SP) resonance of the metalUc particle 112514281 . and 
hmits the SP field enhancement that in the LRA sometimes diverges 12914311 . We employ a 
simple generalization to the LRA, namely the hydrodynamic Drude model (HDM) 12011211 . 
which takes the nonlocal response into account. New in this model, as compared to the LRA, 
are longitudinal waves with sub-nanometer wavelengths, besides the usual transverse waves. 

Incidentally, some HMMs have been found to exhibit strong effective 'nonlocal response', 
even in studies that employ the local -response approximation 1T514T91 . What is meant here is 
that the HMMs, when considered as scatterers, cannot be described in the single-scattering 
Born approximation. This differs from the (material) nonlocal response of HMMs that we con- 
sider here and in Ref. [32], where already the scattering potentials associated with the metal 
constituents are nonlocal. 

We recently showed that material nonlocal response plays an important role on HMMs in the 
limit of vanishing unit-cell size l32l . In particular, the nonlocal response gives rise to a cutoff to 
the hyperbolic dispersion curve, and an associated finite but very large fundamental upper limit 
to the enhanced local optical density of states (LDOS), which in the LRA is known to diverge in 
the limit of vanishingly small unit cells. Now in realized HMMs, it is predominantly the finite 
size of the unit cell that keeps the LDOS finite. So it is an intriguing question, not explored yet 
as far as we know, whether nonlocal response can also have noticeable effects in state-of-the-art 
HMMs. This theoretical paper describes our search for observable nonlocal effects in HMMs 
with unit cells much smaller than an optical wavelength, but much larger than the wavelength 
of hydrodynamic longitudinal pressure waves. 



2. Dispersion relations of hyperbolic metamaterials 

As the HMM we consider a ID sub wavelength dielectric-metal Bragg grating, with a unit cell 
of thickness d, and thicknesses a and b of the dielectric and metal layers, respectively. The 
permittivity of the dielectric layer is ej. The metal is described in the HDM as a free-electron 
plasma with 
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where is the Drude permittivity for the transverse electric fields as in the LRA, while the 
wavevector dependence of the permittivity of the longitudinal electric fields is responsible 
for the nonlocal response. 

We apply the hydrodynamic generalization of the common transfer-matrix method for lay- 
ered systems 02211321 . and obtain the exact dispersion equation of the HMM 
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where for convenience we introduced the dimensionless parameters 
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and where k^^ represents the Bloch wavevector in the direction of the periodicity, k\)^ the wavevec- 
tor along the layers, and ko ^ co/c the free-space wavevector We also introduced the derived 

wavevectors = \Jkled-kf^, kj^ = ^kle'^-kj, and k^^ = y^(a)2 + iyco - a)^)/l5^ - fc^. 

We note that solving the hydrodynamic Drude model requires boundary conditions in ad- 
dition to the usual Maxwell boundary conditions. The dispersion relation Eq. ^ was derived 
using the continuity of the normal component of the free-electron current as the additional 
boundary condition. Details can be found in Refs. Il25ll32ll33l . The exact dispersion relation 
for nonlocal response Eq. (|2) reduces to the exact dispersion relation in the LRA by putting 
the longitudinal permittivity e^(A:, (o) of Eq. ([T]) equal to the transverse permittivity £^(0) of 
Eq. ([T]i, in other words by taking the nonlocal parameter (5 to zero. 

3. Effective nonlocal material parameters 

As stated before, we consider sub-wavelength unit cells much smaller than optical wavelengths, 
but with metal layers much larger than the wavelengths of their longitudinal pressure waves. So 
we focus on the situation where the following two parameters are small, 

M«l, (4) 



where k^ = ^ O)^ + iyco — CO^ / p is the longitudinal wavevector Furthermore, we will only 

consider the frequency range (O < Op, where < and the dispersion curve of the HMM 
could be a hyperbola in the LRA. We now make a first-order Taylor approximation in the 
small parameters of Eq. (|4) to the exact dispersion relation Eq. and obtain the approximate 
dispersion relation 

gioc r gioc ^ 
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The two terms in the brace are two leading correction terms that originate from the nonlocal 
response of the free electrons and from the finite size of the unit cell, respectively. Without 
them, i.e. when neglecting both the finiteness of the unit cells and nonlocal response, we have 
the well-known dispersion relation in the LRA 

gloc 

^i = ^oer-^?i^, (6) 
where e||°'^ and ej"*^ are the effective local permittivities ||9] 

£|| = /d£d +/m£mj £± ~ ~ ~ ~ ~^ (7) 

Obviously, when e||°'^e|°'^ < 0, Eq. ^ describes a hyperbolic dispersion curve. However, we 
are now rather interested in nonlocal response effects in reaUstic HMMs, so we investigate the 



importance of the terms in the brace in Eq. (|5), which are given in terms of 



ghdm ^ and (8a) 

Ala = J2^d^mfdfmd^ i ^0 ' 



T 




12 II 



tdtm td t-r 



(8b) 



From the expression (|8a] i for g^"*™, it follows that the local-response limit is found in the limit 
k^^d °o rather than by taking the size d of the unit cell to zero, see also Ref. 1321 on this point. 
In the following we will keep the first [i.e. the nonlocal) correction term, while neglecting the 
other correction term Ai^. The sole justification for doing so is that the second term turns out 
to be negligible in comparison to the first one for the HMMs that we consider, as we illustrate 
numerically below. So, as an improved approximation to Eq. we keep the correction term 
concerning ej"^™, and then obtain the approximate nonlocal dispersion relation that is central to 
our present study, 

gloc 

kl-klelr-ki-L^, (9) 

in terms of the new effective permittivity in the direction of the periodicity corrected by the 
nonlocal response 

' ' ' (10, 
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The other tensor component ej°'^ of the dielectric tensor has no such nonlocal correction. 

Since l/\k\^d\ ^ 1, we usually have el*"^ ^ ej"^™, and then e"'°'^ « el"*^ is a good approxima- 
tion. However, for large e|°'^ and especially in the extreme case that e^f^ — °°, it may occur that 
e^f^ becomes much larger than e^''™, so that £"'°'^ « ej"^™ by virtue of Eq. (fTOl i. So in this case, 
the nonlocal response has the important role of replacing the infinite el"*^ by the finite e^''™ in 
the dispersion relation. But when does it occur that e|°'^ — !> oo? This immediately follows from 
Eq. (|7J, and occurs when fdE^+fmSd vanishes. In particular, neglecting the Drude damping 7 
of the metal, the frequency 0'°^ for which e|°'^ — > 00 is given by 



Let us also consider another extreme case, namely £"'°'^ — > 00. This is achieved when e|°'^ = 
_ghdm^ which happens at the frequency 
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At this frequency (0^^°'^, the nonlocal response changes the finite e|°'^ into the infinite e"'"*^, 
again a strong nonlocal effect. 

From the above analysis, it follows that nonlocal response is important for HMMs near the 
frequencies (ol°^ and (0^^°'^, as also evidenced numerically in the next section. The larger the 
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difference between col°^ and Ci)j" s , the broader the frequency range with noticeable nonlocal ef- 
fects. As indicated by Eqs. (fTTT i and (fT2l i. the nonlocal resonance frequency (0^^°^^ is blueshifted 
with respect to the local one (ol°^. The relative blueshift {(0"^°'^ — (ol°^)/(ol°^ is proportional 
to Cd, and inversely proportional to the dimensionless parameter k^^a. The latter dependence 
is rather surprising, since k^^ is the longitudinal wavevector of the free electrons in the metal 
layers, whereas a is the thickness of the dielectric layer! We checked, also numerically, that our 
first-order Taylor expansion of the dispersion relation in the small parameter (k^^b)^^ indeed 
gives a relative blueshift ed/(A:J^a), independent of the thickness b of the metal layer. 

4. Effects of nonlocal response on the dispersion curve: numerical analysis 

To numerically illustrate the effects of the nonlocal response on the HMM, we choose a specific 
example of the HMM with a = 6nm, b = 3nm, and = 10. We choose the metal to be Au, and 
describe it by only its free-electron response, with parameters hcOp = 8.812eV, hy = .07 SleV , 
andvF= 1.39 X lO^m/s. 




Fig. 1. (a) kod and (b) real part of e|'|°'=, and (c) real parts of e|°'^ and e^"'™, of the 

HMM. The unit cell of the HMM shown in the inset of (a) has a = 6nm, b = 3nm, = 10, 
and the yellow metal layer is Au. 

Fig.llfa) shows the value of k^d and 1 /|A:[;;/7|. Clearly, in the frequency range O.lcOp < (0 < 
0.6(0p, both kod and 1 satisfy Eq. (|4|. This indicates that we are in the regime where the 

analytical results of the above section are valid. Fig. [Tib) displays the real part of e'"'^ — e"'"'^, 

and Fig. [He) the real parts e^f^ and e"'"*^. The two dashed vertical lines in Fig. [Ttc) mark the 
positions of (ol"^ and O^g"'^ as calculated with the approximate but accurate Eqs. (flU and ( fT2b . 
respectively. With the damping loss of the metal, e^'' and e"'°'^ show large but finite values 
near (ol°^ and Co"^°'^, respectively. Around to'°^ and C0^^°'^, el"*^ and e"'"*^ show a huge difference, 
indicating the importance of nonlocal response, consistent with the theoretical prediction in 
the previous section. Far from (ol°^ and (0^°^, the dielectric functions e|°'^ and e"'°'^ are nearly 
equal, so that away from these resonances the nonlocal response is only a small perturbation. 

Fig. 121 a) and (b) illustrate the dispersion curves at O) = 0.1 (Op and co = 0.6(0p, respectively. 
These frequencies are far away from the resonances (ol°^ and (0"^°^. The small damping loss is 
neglected for a clearer illustration of the dispersion curve. Several observations can be extracted 
from Fig. ID Firstly, the exact local and nonlocal dispersion curves agree well with those based 
on the effective material parameters, which confirms the validity of the approximate theoretical 
results of Sec.|3] Secondly, the local and nonlocal curves only show a slight difference, agreeing 
well with Fig. [Tic) where we see e"'"*^ « el*"^ both for ft) = O.lcOp and 0.6ft)p. 

Fig l3a) and (b) show the on-resonance dispersion curves, namely for ft) = (O^^ = QAlcOp = 
and for ft) = 0)"^°'^ = 0.47 ft)^, respectively. Loss is again neglected, as in Fig. [2] At both 




Fig. 2. Dispersion curves of the HMM at (a) CO = O.lcOp and (b) ft) = 0.6a)p. The HMM 
is as that in Fig. [T] except that the loss is neglected. Red solid curves for the exact local 
dispersion, blue solid curves for the exact nonlocal dispersion, red dashed curves for the 
approximated local dispersion of Eq. (|6), and blue dashed curves for the approximated 
nonlocal dispersion of Eq. (|9}. 
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Fig. 3. Dispersion curves of the HMM at (a) ft) = co}^ = 0.41 ftjp and (b) m = cOj^'f^ = 
0.47 ftjp. The HMM is as that in Fig.[T]except that the loss is neglected. Red solid curves de- 
note the exact local dispersion, blue solid curves the exact nonlocal dispersion; red dashed 
curves for the approximated local dispersion curves of Eq. and blue dashed curves for 
the approximated dispersion curves of Eq. (|9). 



frequencies, the nonlocal response modifies the dispersion curve noticeably. In particular, at 
CO = 0.41a)p, the local dispersion curve consists of two nearly flat lines of fc^ = iko^Je^. 
With the nonlocal response, however, the dispersion curve becomes a closed ellipse! This re- 
markable difference can be understood from the difference between e^f^ and e"'*"^. In particular, 
gioc diverges while e"'"*^ = 41 stays finite in the lossless case. When including the loss, the ef- 
fective parameters change to e^f^ = —200 + 500/ and e"'°'^ = 41 +2.5/. The loss would only 
slightly modify the dispersion curves of Fig. |3ja). In Fig. |3jb) we display dispersion relations 
at o = 0.47 (Bp. Here the local dispersion curve is a hyperbola. With the nonlocal response, how- 
ever, the dispersion curve becomes nearly flat lines, which can be attributed to the extremely 
large value of e"'°'^ at this frequency CO^gf^. 

5. Effects of nonlocal response on a hyperbolic metamaterial lens 

A HMM slab can operate as a superlens with subwavelength resolution, since the hyperbolic 
dispersion curve supports propagating waves with arbitrarily high wavevectors, which can 
transfer the evanescent information of the object. In the LRA, it is known that the HMM with its 
flat dispersion curve at a)'°s [recall Fig. ^a)] is especially favorable for subwavelength imag- 
ing i8J. The reason is that the flat dispersion curve with the constant k± ensures that all plane- 



wave components experience the same phase changes after transmission through the HMM 
slab, at least if reflections can be neglected. Actually, the reflections can be suppressed by ap- 
propriate choice of the thickness / of the HMM lens, and even be made to vanish by choosing 
k±l = iiTZ, where n is a positive integer In theory this could lead to a perfect image at x = 0. 

Let us now investigate how nonlocal response may influence the subwavelength imaging 
characteristics of the HMM lens. As demonstrated in the above sections, the nonlocal response 
sets the infinite e|°'^ ^ 0° to a finite value and accordingly destroys the desired flat dispersion 
curve at the frequency col"^ where one would normally choose to operate for perfect imaging, 
based on the LRA. Nevertheless, at a blueshifted frequency (o"^°^, the local effective dielectric 
function e^f"^ may be finite, but the hydrodynamic Drude model predicts instead that £"'°'^ di- 
verges, with the concomitant flat dispersion curve suitable for subwavelength imaging. Thus 
nonlocal response is expected to strongly affect the performance of HMM superlenses for fre- 
quencies (0 around (ol°^ and C0"^°^. 
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Fig. 4. Transmitted electric-field intensity distribution for a line dipole source J = S{x + 
I +Xs)8{y)y positioned to the left of a HMM slab that has its left interface alx = ~l and 
right interface at x = 0. For co = OAlcOp ~ cojes ' panel (al) shows the intensity for local 
response, (bl) for nonlocal response, and (cl) the intensity along y with x = 0. Panels (a2), 
(b2), and (c2) are the analogous graphs for co = OAlcOp ~ co"g^- In panels (cl) and (c2), 
the green dashed curves are for the case without the HMM slab, the red solid curves for the 
HMM with local response, and the blue solid curves for the HMM with nonlocal response. 
The HMM unit cell is as that in Fig.[TJ x^, = Inm and / = 36d. 



As an example, we consider a HMM slab with I = 36d, i.e. composed of 36 unit cells, in 



a free-space background with the two boundaries at x = —I and x = 0. The unit cell is as in 
Fig. [T] and is arranged in a symmetric sandwich structure with the metal layer at the center. 
A line dipole source is positioned to the left of the HMM slab, with jc-coordinate — / — Xs and 
y-coordinate 0, and is represented by J = 8{x + l -\- x^)5{y)y. We choose the distance to the 
HMM slab to be x^ = Inm. 

Fig-Sal) and (bl) demonstrate the transmitted electric-field intensity at © 0)'°^ = 0.41a)p 
for local and nonlocal response, respectively. Metal loss is taken into account. In the local case, 
we have the known the dispersion of Fig. [Sja) featuring flat lines. Additionally, the reflection 
is nearly zero since kj k, Ak, so all wave components experience the same phase change. 
Accordingly, after propagation through the HMM from x = — / to x = 0, a subwavelength image 
is formed near x = 0, as seen in Fig. Sfal). By contrast, in the nonlocal case the dispersion 
curve turns into an ellipse. This closed dispersion curve sets an upper wave vector cutoff for 
the evanescent waves, and the wave components below the cutoff experience different phase 
changes. Accordingly, the quality of the image becomes worse, see Fig.|4ta2). 

In Fig.Scl), we depict the electric-field intensity again at (o\°l as a function of y at the HMM 
boundary x = 0, and also for the case without the HMM slab. It is seen that the electric-field 
intensity is nearly vanishing in the absence of the HMM slab, since the evanescent components 
vanish after traveling a distance of /, which is of the order of one free-space wavelength. With 
the HMM slab in place, the electric-field intensity in the local case shows a subwavelength 
image that peaks at y = 0, with a full width at half maximum (FWHM) of only 0.47 Ap. With 
nonlocal response, however, the electric-field intensity distribution for (ol°^ becomes flatter, 
with a double rather than a single peak, with peak intensities at y = ±0.55Ap. 

Let us now turn to the other resonance frequency, namely o = fi)|"'°'^ of Eq. (fT2l i. Fig.Ua2) 
and (b2) show the transmitted electric-field intensity for the local and nonlocal cases, respec- 
tively, at o = ©res"^ = 0.47 Op. At this frequency, the local dispersion curve is a hyperbola, 
while the nonlocal dispersion shows two flat lines. In the nonlocal case, we have kj k, A.Ik in- 
dicating that reflections do exist and the wave components experience different phase changes. 
However, compared to the local case, the transmission coefficients for the different wave com- 
ponents in the nonlocal case vary more smoothly as a function of k||, thanks to the flat disper- 
sion curve with the constant A:^. As a result, the hydrodynamic Drude model predicts a better 
focusing performance of the HMM slab at e) = 0.47 (Bp = (0^1°^ than does the local-response 
theory, compare Fig. 21 a2) and (b2). In Fig. |4jc2), the electric-field intensity along y at x = 
is shown. In the local case, the electric-field intensity shows several peaks with the strongest 
two at y = ±0.33Ap. In the nonlocal case, the electric-field intensity is peaked at y = with a 
subwavelength FWHM of only 0.48Ap. 

6. Detecting nonlocal response by near-field measurement 

Fig.|4]indicates the possibility of detecting the nonlocal response experimentally by measuring 
the transmitted near-field distribution at the surface of a HMM superlens, which in our setup 
would be its right interface at x = 0. The images in Fig. |4|cl,c2) correspond to hypothetical 
measurements with infinitely small detectors. In experiments, the measured near-field signal 
will rather be an area-averaged electric-field intensity 



where D is the finite detection area of the detector Let us now assume that we have a near- 
field detector with detection area in the yz-plane, with a square shape of size 40 nm x 40 nm, 
touching the HMM interface at x = 0. For the same Ught source interacting with the HMM 
slab as in Fig. ID we depict in Fig. |5]the calculated /av as a function of the y-coordinate of 




(13) 



D 



the center of the detector. It is seen that local and nonlocal response models give significantly 
different predictions for the measured signal, and thus the differences between the two models 
survive detection area averaging. The single peaks at OAlcOp for local response and at QAlcOp 
for nonlocal response are naturally broader than in Fig. 5Jc l,c2), also due to the area averaging. 
Interestingly, as the frequency increases from 0.41 fflp to 0.47 fflp, the distribution of /av becomes 
broader in the LRA, but narrower in the hydrodynamic Drude model. 
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Fig. 5. Calculated detector-area-averaged near-field intensity 4v of the light emitted by the 
same source as in Fig. |4] and transmitted through the same HMM slab. The detector is 
square-shaped with size 40nm x 40nm, touching the HMM right interface at x = 0. It scans 
along the y-direction, and the y-coordinate on the horizontal axes in panels (a) for the LRA 
and (b) for the HDM corresponds to that of the center of the detector. 



7. Conclusions 

We investigated the effects of the hydrodynamic nonlocal response on hyperbolic metamate- 
rials with periodicity in the subwavelength regime of the transverse optical wave, but much 
larger than the wavelength of the longitudinal hydrodynamic pressure waves. It is found that 
the nonlocal response corrects the effective permittivity tensor element e|°'^ in the periodicity 
direction to the new form e"'°'^ of Eq. ( fTOl i. Around the frequencies a^"^ and (0^^°'^ correspond- 
ing to el"*^ — > oo and e"'°'^ — !> respectively, e|°'^ and e"'"*^ show noticeable differences, even 
leading to completely different dispersion curves with and without the nonlocal response. 

We find that nonlocal response blueshifts the resonance frequency of HMMs from (o\°l 
[Eq. (fTTT il to (o"l°^ [Eq. (fT2bl. The relative blueshift has an interesting simple form independent 
of the thickness of the metal layers. Similar nonlocal blueshifts for single nanoplasmonic par- 
ticles have been predicted before, and significant blueshifts have also been measured II26II28I : 
how much of these can be attributed to hydrodynamic effects is a hot topic 1 26, 281 1341 . 

Furthermore, we predict that nonlocal response shows its mark in the performance of a finite 
HMM slab as a focusing lens: when increasing the operating frequency from (ol°^ to (0^1°'^, for 
local response the near-field distribution of the transmitted light shows an image that gets out 
of focus, whereas the focus would improve instead according to the nonlocal-response theory. 
We propose to test the blueshift and these contrary predictions experimentally, as a clear and 
interesting test whether nonlocal response can be observed in hyperbolic metamaterials. 
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